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Abstract 

Percolation models with multiple percolating clusters have attracted much attention in recent 
years. Here we use Monte Carlo simulations to study bond percolation on L\ x L2 planar random 
lattices, duals of random lattices, and square lattices with free and periodic boundary conditions, in 
vertical and horizontal directions, respectively, and with various aspect ratio L\jLi- We calculate 
the probability for the appearance of n percolating clusters, W n , the percolating probabilities, P, 
the average fraction of lattice bonds (sites) in the percolating clusters, < c b > n (< c s > n ), and the 
probability distribution function for the fraction c of lattice bonds (sites), in percolating clusters 
of subgraphs with n percolating clusters, f n (c b ) (fn{c s ))- Using a small number of nonuniversal 
metric factors, we find that W n , P, < c b > n (< c s > n ), and f n (c b ) (f n {c s )) for random lattices, 
duals of random lattices, and square lattices have the same universal finite-size scaling functions. 
We also find that nonuniversal metric factors are independent of boundary conditions and aspect 
ratios. 

PACS numbers: 05.50.+q, 64.60.Ak, 75.10.-b 
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I. INTRODUCTION 



Percolation is related to many interesting scientific phenomena [JT]. In recent years per- 
colation problems with multiple percolating clusters have attracted much attention [2-19]. 
Most of the simulational studies of such problems have been restricted to percolation on 



lattices [20|. However, many physical systems with multiple percolating clusters such as 
Carbino disks used in the study of quantum Hall effects @], or oil fields confronted with 
drilling problems, do not have underlined regular lattice structures. Thus, it is of interest 
to know the relationship between the quantities for percolation on regular lattices and the 
quantities for percolation not on regular lattices, such as random lattices. In the present 
paper, we use Monte Carlo simulations to study bond percolation on L\ x L2 planar random 
lattices, duals of random lattices, and square lattices with free and periodic boundary condi- 
tions in vertical and horizontal directions, respectively, and with various aspect ratio L\jLi. 
We calculate the probability for the appearance of n percolating clusters, W n , the percolating 
probabilities, P, the average fraction of lattice bonds (sites) in percolating clusters, < c b >„ 
(< c s > n ), and the probability distribution function for fraction c of lattice bonds (sites), 
in percolating clusters of subgraphs with n percolating clusters, f n (c b ) (f n (c s )). Using a 
small number of nonuniversal metric factors, we find that W n , P, < c b >„ (< c s > n ), and 
fn{c b ) (f n {c s )) for random lattices, duals of random lattices, and square lattices, have the 
same universal finite-size scaling functions. We also find that nonuniversal metric factors are 
independent of boundary conditions and aspect ratios. Furthermore, this study is related 
to recent developments in the universality and scaling of critical phenomena. 

Universality and scaling are two important concepts in modern theory of critical phenom- 
ena [21], |23| , and percolation models are an ideal system for studying critical phenomena 
IJ. Thus, universality and scaling have been actively studied in recent decades, especially 
for percolation models pi |. In 1992, Langlands et. al. |25] proposed that for bond and site 
percolation models on square (sq), planar triangular (pt), and honeycomb (he) lattices, the 
critical existence probability (also called crossing probability or spanning probability) is a 
universal quantity, when aspect ratios of sq, pt, and he lattices have relative ratios, 1 : \/3 : 
y / 3/2. From 1995~1996, Hu, Lin and Chen (HLC) H |26| calculated existence probability, 
E p , percolation probability, P, and probability for the appearance of n percolating clus- 
ters, W n , of bond and site percolation models, on sq, he, and pt lattices with aspect ratios, 



2 



1 : v^3 : v / 3/2, and showed that all their scaled data fall on the same universal scaling 
functions, by selecting a very small numbers of nonuniversal metric factors, and maintain- 
ing similar nonuniversal metric factors under two boundary conditions, free and periodic 
boundary conditions. By using renormalization group theory, Hovi and Aharony in 1996 
p7| also pointed out that scaling functions, for the spanning probability are universal at 
the fixed point for every system with the same dimensionality, spanning rule, aspect ratio 
and boundary conditions. Okabe and Kikuchi in 1996 fl28|, extended the work of HLC to a 



two-dimensional Ising model on planar regular lattices. In 1997, Hu and Wang [ 1 1 found 
that the lattice percolation and continuum percolation of hard and soft disks, have the same 
universal scaling functions for W n . Using the connection between an Ising model and a 
bond-correlated percolation model P5fl , Tomita, Okabe, and Hu in 1999 |I7|, calculated the 
probability for the appearance of n percolating clusters, W n , the percolating probabilities, 
P, the average fraction of lattice sites in percolating clusters, < c > n , and the probability 
distribution function for the fraction c of lattice sites in percolating clusters of subgraphs 
with n percolating clusters, f n (c) , for bond-correlated percolation model on sq, he, and pt 
lattices, with aspect ratios of 1 : v^3 : y/3/2. Using a small number of nonuniversal metric 
factors, they found that W n , P, < c > n , and f n {c) for sq, he, and pt lattices have the same 
universal finite-size scaling functions. 

However, the studies mentioned above are mostly focused on regular lattices, with fixed 
coordination numbers |[20|| . In 1999 Hsu and Huang (HH) |5(J determined the percolation 
thresholds and critical exponents, and demonstrated explicitly that the ideas of universal 
critical exponents and universal scaling function with nonuniversal metric factors can be 
extended to bond percolation on L x L periodic planar random lattices, duals of random 
lattices, and square lattices, for both existence and percolating probabilities, and mean 
cluster size. This paper will study bond percolation on L\ x L 2 planar random lattices, 
duals of random lattices, and square lattices, in more detail, and consider the case that the 
lattices have free and periodic boundary conditions, in vertical and horizontal directions, 
respectively, as in 0. Percolating probability is defined by the ratio of the number of 



bonds in the percolating clusters to the total number of bonds in ||30|| . Here, we consider 
two different definitions of the percolating probability, in terms of bonds and in terms of 
sites; the latter was also used in [Q and JH]]. We calculate the probability W n for the 
appearance of n percolating clusters, the percolating probability P, the average fraction 
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of lattice bonds (sites) in percolating clusters, < c b > n (< c s > n ), and the probability 
distribution function for the fraction c of lattice bonds (sites) in percolating clusters of 
subgraphs with n percolating clusters, f n (c b ) (/ n (c s )), for various values of aspect ratios 
Lx/L 2 , and finally check the universal finite-size scaling behaviors for these quantities. In 
pOfl , HH used two nonuniversal metric factors, D 2 and D%, to fix universal finite-size scaling 
functions, for percolating probability in terms of bonds. In the present paper, we calculated 
two nonuniversal metric factors of percolating probability, in terms of sites and obtained 
previous known values of nonuniversal metric factors determined by HH, to check whether 
we could have universal scaling functions for W n , P, < c b > n (< c s > n ), and f n (c b ) (f n (c s j), 
of bond percolation on random lattices, duals of random lattices, and square lattices. 



Dirichlet and Voronoi [[31]] first used the concept of random lattices in condensed matter 
theory and Christ, Friedberg, and Lee (CFL) |32|] used another type of random lattices to 
formulate quantum field theory. Here, we adopt the CFL algorithm and give a brief review 
of the construction of planar random lattices and their duals. First, we randomly generate 
N sites in the Li x L 2 rectangular domain with periodic boundary conditions. Next, we 
choose three nearby sites arbitrarily, and draw a circle to go through the three sites. If there 
are no lattice sites inside the circle, the three sites are connected by links to form a triangle. 
A planar random lattice is constructed by repeating the process until all sites are connected 
by links. The whole rectangular domain is divided into 2N non-overlapping triangles, whose 
vertices are sites of the random lattice, and circle centers with triangles are the sites of dual 
lattice. Thus, there is a one to one correspondence between triangles and dual lattice sites. 
Because a link of the random lattice is shared by two triangles, the two corresponding dual 
lattice sites are connected by one dual link. There is a one to one correspondence between 
links and dual links. The whole rectangular domain is partitioned into N non-overlapping 
planar convex polyhedra, which are formed by dual links and the vertices of iV polyhedra 
are sites for dual lattice. There is also a one to one correspondence between the lattice sites 
and polyhedra on dual lattice. Examples of a planar random lattice with its dual, under 
periodic boundary conditions, in both vertical and horizontal directions, are shown in Fig. 
1(a). 

This paper is organized as follows: In Sec. II, we present the simulational results for 
W n , P, < c b > n (< c s >„), and f n (c b ) (f n (c s )) for bond percolation, on L\ x L 2 random 
lattices, duals of random lattices, and square lattices, under free and periodic boundary 
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conditions in vertical and horizontal directions with Li/L 2 = 4. The boundary bonds 
which cross the rectangular domain in vertical direction on random lattices, due to periodic 
boundary conditions, are eliminated because of free boundary conditions in vertical direction 
considered in this paper. We adopt the method of HH |30| , to find percolating clusters. Only 
the first kind of percolating cluster paths, without boundary bonds in the vertical direction 
(the clusters extend from top to bottom), should be identified, and an example of this 
is shown in Fig. 1(b). In Sec. Ill, we use finite-size scaling theory to check the scaling 
behaviors of various quantities, and to show that such quantities have universal finite-size 
scaling functions for regular lattices and random lattices. A summary is provided in Sec. 
IV. 



II. W n (L ll L 2 ,p), f n (c), and < c>„ 

We see the bond percolation on a lattice G, with linear dimensions L\ and L 2 in horizontal 
and vertical directions, respectively; the probability for the appearance of n top-to-bottom 
percolating clusters, W n (Li, L 2 ,p), is defined by || 

w n {UM,v) = £ P biG ' n) ^ -p) E - b{G ' n) - (i) 

G' n CG 

Here, the percolating cluster is defined as a cluster extending from top to bottom in G, G' n 
denotes a percolating subgraph with n percolating clusters, b(G' n ) is the number of occupied 
bonds in G' n , and E is the total number of links in G. The existence probability E p can be 
expressed obviously as 

oo 

E p = Y,W n , (2) 

n=l 

with Wo = 1 — E p . 

To obtain more detailed information about the contents of the percolating cluster, fol- 
lowing Tomita et. al. [T7|, we decompose W n as 



»i 

W n = I f n {c)dc, (3) 



o 



where n — 1, . . . oo, c denotes the fraction of lattice bonds (sites) in percolating clusters, and 
f n (c) is the probability distribution function of c in subgraphs with n percolating clusters. 



5 



The probability distribution function of c in all subgraphs is the overall summations of f n (c), 
i.e., 

oo 

f(c) = J2fn(c). (4) 

n=l 

In terms of f n (c), the average fraction of lattice bonds (sites) in subgraphs with n percolating 
clusters can be expressed as 

< c > n = / cf n (c)dc, (5) 
Jo 

where n — 1, . . . , oo, and percolating probability P can be written as 

x poo 

< c >= ^2 < c >n= / cf(c)dc = P. (6) 

To generate subgraphs, we use the random bond occupation process with equal occupation 
probability for each link. The simulations are performed on 128 x 32, 256 x 64, and 512 x 128 
planar random (pran) lattices, and their duals (dpran), with free and periodic boundary 
conditions in the vertical and horizontal directions, respectively. To compare the results with 
regular lattices, we also perform simulations on square (sq) lattices of the same sizes. On 
each lattice, we take 60 occupation probabilities around the critical percolation threshold for 
every 0.002 increment, and use the random bond occupation process to generate 10 5 ~ 10 6 
configurations for each occupied probability, p. We calculate W n (Li, L 2 , p) , < c b > n , and 
< c s > n , where c b denotes the fraction of bonds in percolating clusters, and c s denotes the 
fraction of sites in percolating clusters; and the results are shown in Figs. 2 and 3. The 
calculated results of the percolating probabilities in terms of bonds, P b , and in terms of sites, 
P s , are also shown in Fig. 3. We calculate f n (c b ) and f n (c s ) at p = p c and take p c = 0.3333 
for planar random lattices, and p c = 0.6667 for dual lattice p0|| . The results are shown in 
Fig. 4. The differences between bond and site contents in percolating clusters are shown in 
Figs. 3 and 4; here, for the clarity of presentation only, the results for 512 x 128 lattices are 
plotted in the figures. 



III. UNIVERSAL FINITE-SIZE SCALING FUNCTIONS 



The finite-size scaling theory was first formulated by Fisher in 1971 p^j . According to 
the theory, for a physical quantity X, which scales as X(t) ~ t p in a thermodynamic system 
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near critical point t — 0, then the same quantity in a finite system with linear dimension L, 
Xi(t), should obey the general law, 

X L (t) ~ L- p l v F{tl}l"). (7) 

Here, F(x) with x = tL 1 ^ labelled as scaling function, with v as correlation length ex- 
ponent. In 1984, Privman and Fisher considered universal finite-size scaling functions 
and nonuniversal metric factors, and proposed that the singular part of the free energy of a 
critical system can be written as 

f' L {t) ~ L- d Y(DtL^), (8) 

where d is the spatial dimensionality of lattice, Y is a universal scaling function, and D is a 
nonuniversal metric factor. 

At the critical point, p = p c , there also exists a finite-size scaling form for the distribution 
function of X L (t) 

Q(X L (t = 0)) ~ L"' v Y{X L {t = 0) ■ L p/v ). (9) 



In p6| , |30| , three nonuniversal metric factors Di, D 2 , and D 3 were used for regular lattices 
and random lattices, to describe the universal scaling functions of existence probability E p , 
and the percolating probability P, i.e. 

E p (p,L) = F(x), (10) 

with x = Di(p — p^L 1 ^ and 

D 3 P(p,L) = L-^S p (z), (11) 

with z = D 2 {p - p c )L 1/v '. 

Following Hu et. al. 0, we use the evaluated percolation thresholds p c [pOj] , and the 
exact values of critical exponent z/ = 4/3, to plot W n as a function of y = (p — p^L 1 ^, 
for planar random lattices and their duals in Figs. 5(a) and 5(b), respectively. We can see 
from these results, that the scaled data for W n can be described by a single scaling function 
F n (r,y) with r = Li/L 2 , and F n (r,y) for n > 2 as a symmetric function of y. In Fig. 6 
we plot W n (Li, L 2 ,p) as a function of x for bond percolation on 512 x 128 random lattice, 
dual of random lattice, and square lattice, where x = D\{p — p c )L^ v , with D\ taken from 
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Table I and L = (Li x Lq) 1 ! 2 . Fig. 6 shows that the calculated results for each n can be 
well described by a single universal scaling function, U n (x). 

In Hu and Lin's paper M, the scaling functions F n (r,y) were calculated for bond per- 
colation on square lattice for various values of aspect ratios, r. We will examine whether 
the same nonuniversal metric factors, Di, can be extended to different aspect ratios. We 
calculate W n for L\ x L 2 random lattices, dual of random lattices, and square lattices with 
r = L1/L2 = 1, 2, . . . , 6 and determine the universal scaling functions, U n (r,x), where 
x = Di(p — p c )L l l v and where D 1 is taken from Table I. The results for n — 1 and 2 are 
shown in Figs. 7(a) and 7(b), respectively. We can see that the scaled data for each r can 
be described by a single universal scaling function very well. The results of U n ( t ', x) as a 
function of r for n — 0, 1, . . . , 4 at the critical point p = p c , are presented in Fig. 8(a) which 
show that the three lattices provide similar results. We also calculate the average number 
of percolating clusters C(r,x), defined by 

00 

C(r,x) = J2 U n(r,x)n. (12) 

n=l 

C(r, 0) for random lattice, dual of random lattice, and square lattice as a function of r are 
shown in Fig. 8(b). Fig. 8(b) shows that C(r, 0) increases linearly with an increasingly 
large, r and that different lattices have the same slope of approximately 0.43. 

Tomita et. al. [T7| had obtained universal finite-size scaling functions of < c > n , < c >, 



/ n (c), and f(c) for a bond-correlated percolation model, corresponding to an Ising model on 
planar regular lattices. It is of interest to extend such study to bond random percolation on 
random lattices. From Eqs.(5), (6), and (11), the universal scaling function of < c > and 
< c > n can be expressed as 

D 3 <c(jp,L)>=L-V v G(z), (13) 

and 

D 3 <c(p,L)> n =L-^G n (z), (14) 

with z = D 2 (p — Pc)L l ' v . At p = pc, the universal scaling function of /(c) and / n (c) are 
expressed as 

D^f(c)=L^H(z>), (15) 
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and 



Dz 1 f n (c)=L^H n (z'), (16) 

with z' = DscL 13 ^. To check the finite-size scaling and universality of these quantities, we 
use simulation results for 256 x 64 and 512 x 128 square lattices, planar random lattices and 



their duals. In |3(| , the percolating probability P is defined in terms of the bond number in 
percolating clusters, and the nonuniversal metric factors D 2 = D\ and D% = D b are used. 
To evaluate factors D s 2 and £)|, we adopt the same procedure as in pD| , plotting P b /L~ )3 / U 



and P s /L $' v as functions of y = (p — p^)L l ' v as shown in Fig. 9. All of the nonuniversal 
metric factors for different types of lattices used in this paper are listed in Table I. 

We plot D 3 P/L-M V and D 3 < c > n /L~< 3 I V as a function of z = D 2 (p - p c )L 1/u , in 
Figs. 10(a) and 10(b) for bond content and site content, respectively, with D 2 and D 3 taken 
from Table I. At p = p c , the scaled data D^ 1 /(c) / "LPI V , and D 3 l f n (c) / L^l v as functions of 
are presented in Figs. 11(a) and 11(b), respectively, for the bond content and 
site content. Figs. 10 and 11 show that the bond percolation processes on square lattices, 
random lattices and their duals have universal finite-size scaling functions. 

IV. SUMMARY AND DISCUSSION 

Having used nonuniversal metric factors from Table I in this paper, we have found that 
the universal finite-size scaling functions for W n (Figs. 6 and 7), < c b > n and P b (Fig. 
10(a)), < c s > n and P s (Fig. 10(b)), f n (c b ) (Fig. 11(a)), and f n (c s ) (Fig. 11(b)). Fig. 7 
includes results for different aspect ratios r, i.e. 6 > r > 1. The values of nonuniversal 
metric factors, D\, D 2 , and D 3 of Table I are consistent with the corresponding values of 



Ref. |30[, where the boundary conditions are different from the boundary conditions of the 
present paper. These results suggest that in random lattices, the nonuniversal metric factors 
are also independent of the boundary conditions and aspect ratios, as in the case of regular 
lattices [p6[| . Please also note that D 2 and _D| in Table I are consistent within numerical 
errors. 

Many interesting problems are related to the properties of multiple percolating clusters. 
It is of interest to extend the study of the present paper to higher spatial dimensions. In 
particular, the further study of multiple percolating clusters in three dimensions could be 
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related to an oil drilling problem. 
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Table I: The values of metric factors D 1: D\, D\, D s 2 and _D|, for square lattices, random 
lattices and their duals, with free and periodic boundary conditions in vertical and horizontal 
directions, respectively. 



Lattices 


Square 


Planar random 


Dual of planar random 


D 1 


1 


1.166 ±0.020 


1.177 ±0.016 


D\ 


1 


1.164 ±0.014 


1.176 ±0.015 


D\ 


1 


1.512 ±0.008 


0.778 ± 0.002 




1 


1.186 ±0.012 


1.180 ±0.014 




1 


1.062 ±0.001 


1.005 ±0.002 
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FIGURES 



FIG. 1. Examples of (a) L\ x L 2 = 8 x 4 planar random lattice (solid lines), with its dual 
(dashed lines) on L\ x L 2 = 8 x 4 rectangular area, with periodic boundary conditions, and 
(b) a first kind of percolating cluster path, without boundary bonds (bold solid lines) on 
random lattice. 

FIG. 2. W n (Li, L 2 ,p) on square lattices, planar random lattices and their duals of size 
128 x 32, 256 x 64, and 512 x 128, which are represented by dotted, dashed, and solid lines, 
respectively. 

FIG. 3. (a) P b and < c b > n , and (b) P s and < c s > n on square lattice, planar random 
lattice and its dual of same size 512 x 128. 

FIG. 4. At p — p c , (a) f(c b ) and f n (c b ), and (b) f(c s ) and f n (c s ) on square lattice, planar 
random lattice and its dual of same size 512 x 128. 

FIG. 5. The scaled results of F n (r,y) = W n (LiL 2) p) as a function of y = (p — p<?)V~l v for 
(a) planar random lattices and (b) their duals of size 128 x 32, 256 x 64, and 512 x 128. The 
monotonic decreasing function is for F (r,y). The S shaped curve is for Fi(r,y). The bell 
shaped curves from top to bottom are for F n (r,y), with n as 2, 3, and 4, respectively. 

FIG. 6. The scaled results of U n (x) = W n (L 1: L 2 ,p), as a function of x — Di(p — p c )L l l u for 
square lattice (solid curves), planar random lattice (dashed curves), and the dual of planar 
random lattice (dotted curves) of same size 512 x 128. 

FIG. 7. U n (r,x) for square lattice (solid curves), planar random lattice (dashed curves), and 
the dual of planar random lattice (dotted curves), (a) n — 1, the intersection of the curves 
on x = axis, from up to down are for r = L ± /L 2 = 1, 2, 3, . . . , 6. (b) The width of curves, 
from small to large are for r = L\jL 2 = 1, 2, 3, . . . ,6. 

FIG. 8. (a) U n (r, 0) as a function of r = Li/L 2 , for a number of percolating clusters (npc) 
run from to 4, and (b) C(r, 0) as a function of r = Li/L 2 , with slope of the fitting line 
0.43. Square lattice (□), planar random lattice (A) and the dual of planar random lattice 
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(x), all have horizontal periodic boundary conditions. 



FIG. 9. The scaled results of P(p, L) /L~^ u , in terms of bonds (P = P b ) and sites (P = P s ), 
for different types of lattice size 256 x 64 and 512 x 128, as a function of y = (p — p c )L l l u . 

FIG. 10. The scaled results of G n (z,L) = D 3 < c(p,L) > n /L^ 13 ^, as a function of z = 
D 2 (p — p c )L 1 / iy , for square lattices, planar random lattices and their duals of same size 
256 x 64, and 512 x 128. (a) c = c b , D 2 = D b 2 , D 3 = D b 3 , and (b) c = c s , D 2 = D s 2 , D 3 = D s 3 . 

FIG. 11. The scaled results of H n (z',L) = D 3 l f n (c) / , as a function of z' = D 3 (p - 
PcjL 13 ^ for square lattices, planar random lattices, their duals of same size 256 x 64, and 
512 x 128, with fitting results (solid curves); (a) c = c b , D 3 = D^, and (b) c = c s , D 3 = L>|. 
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